Frequency and damping of the Scissors Mode of a Fermi gas 
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We calculate the frequency and damping of the scissors mode in a classical gas as a function of 
temperature and coupling strength. Our results show good agreement with the main features ob- 
served in recent measurements of the scissors mode in an ultracold gas of ®Li atoms. The comparison 
between theory and experiment involves no fitting parameters and thus allows an identification of 
non-classical effects at and near the unitarity limit. 
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INTRODUCTION 

By tuning the atom-atom interaction in a gas of 
fermions it is possible experimentally to investigate the 
crossover from a Bose-Einstein condensate (BEG) of 
molecules to a BCS superfluid. In the crossover re- 
gion the system is strongly interacting in the sense that 
kF\a\ 3> 1 (unitarity limit), where a is the s-wave scatter- 
ing length for the atom-atom interaction and k^p = 2>n^n, 
with n being the density of the gas. A theoretical descrip- 
tion of this region is challenging, in particular at non-zero 
temperature 

GoUective mode experiments have given a wealth of in- 
sight into the properties of atomic gases, since they often 
provide more detailed information than e.g. thermody- 
namic measurements 0] ■ In a recent experiment [4| the 
scissors mode excitation in an elliptical trap was used to 
characterize the transition between hydrodynamic and 
coUisionless behavior as a function of temperature and 
scattering length. Since both the frequency and attenua- 
tion is affected by the atom-atom interaction, such mea- 
surements can give important information on the prop- 
erties of a fermion system at and near the unitarity limit. 

In this Brief Report we calculate the frequency and 
attenuation of the scissors mode investigated in Ref. 
assuming the temperature to be sufficiently high that the 
gas can be treated as being classical. Since the calcu- 
lation involves no fitting parameters, it can be used to 
identify non-classical features of the observed frequency 
shift and damping as functions of temperature and/or 
interaction strength. 

KINETIC THEORY FOR A CLASSICAL GAS 

Gonsider a two-component (for brevity denoted by 
"spin" with the two values a =]" , |) Fermi gas of atoms 
with mass m in its normal phase trapped in a potential 
Vir) = m(a;^x^ -f ^yV"^ + w^z^)/2. At high tempera- 
tures, in the classical regime, the dynamics is described 
by a semi-classical distribution function /(r, p, t), which 
satisfies the Boltzmann equation. For the scissors mode 
studied experimentally in [J] , the two components of the 



gas move together, and we need only consider one distri- 
bution function / = /| = fi- 

To calculate the frequency and damping of the scissors 
mode, we linearize the Boltzmann equation in terms of a 
small deviation 6f = f — from the equilibrium distri- 
bution f{Y, p) by writing 5f{v, p, t) = /°(r, p)$(r, p, t). 
The linearized Boltzmann equation becomes 

where f = v = p/m and p = —dVjdv. The collision inte- 
gral / is given in Ref. Q with an interaction described by 
s-wave scattering involving only particles with opposite 
spin. The cross section a is 

_ ^■^"^ 

" 1 + {vrajhY ' ^ ' 

where Pr is the relative momentum of the scattering par- 
ticles. The unitarity limit is defined by \a\ — > oo. 

In the hydrodynamic limit, a scissors mode in the xy 
plane is characterized by a velocity field v cx V(a;?/) or 

v-6(y,a:,0), (3) 

where 6 is a constant. To describe this mode we choose 
an ansatz for (5/ of the form 

^(r, p, t) = {cxxy 4- cix'py 4- c^iyy.j, -f C4P2;Py)e"*"*, (4) 

where the Ci are constants and uj is the mode frequency. 
We insert this ansatz into the linearized Boltzmann equa- 
tion ID) and take moments by multiplying by any of the 
terms xy, xpy,ypx and PxPy appearing in $ and subse- 
quently integrating over both r and p. The result is a 
homogeneous set of four coupled equations for the coeffi- 
cients ci, . . . , C4. The frequencies of the collective modes 
are determined by the roots of the corresponding deter- 
minant which yields the equation 

^(c.^-c.D + (-^--.^)(-^-c.,^.)=0. (5) 

Here 

= JujI+ ujI (6) 



2 



is the mode frequency in the hydrodynamic hmit and 

Wcl ^UJx+^^y and UJc2 = \(^x - t^yl (7) 

" The 



the mode frequencies in the coUisionless Hmit 
viscous relaxation rate in ([3]) is given by 



1_ _ / d^rd^ppxPyI[pxPy] 
T J d?'r(Ppp'^p'^f'^ 

In the classical limit, one obtains 



1 N muj^ A-Ko^ 

-,3 — 



dx 



1 + {T/Ta)x^ 



(8) 



(9) 



with uj — ujxUjyUJz, while the characteristic temperature 
Ta is defined by kTa = fi^/ma^. The integral in ^ comes 
from averaging over momentum and space the cross sec- 
tion multiplied by an appropriate weight function. In 
the unitarity limit, (O gives r'^ = iNfi^Lj^ /[15TT{kT)'^]. 
Using © and © we can now calculate the frequency 
and damping rate F ~ — Imcj of the scissors mode as a 
function of temperature and scattering length. In the 
hydrodynamic limit r ^ 0, we obtain from ([5|) 



UJ ~ LUh — IT- 



In the coUisionless limit r ^ cx3, Eq. (O yields 



I 
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(10) 



(11) 



with i = 1,2. This result, F = 1 /4t, for the damping rate 
in the coUisionless limit, obtained by taking moments of 
the kinetic equation, is in fact exact. This can be seen by 
treating the collision integral as a perturbation, following 
the approach used in Ref. Q for the coUisional relaxation 
of trapped bosons above their condensation temperature. 



COMPARISON WITH EXPERIMENT 

We now compare our calculated frequency and damp- 
ing rate of the scissors mode with the recent experiment 
on trapped ^Li atoms [J]. In Fig. [U we plot the observed 
frequency and damping rate as a function of temperature 
at unitarity, l/kpa = 0. The scissors mode is excited in 
the xy plane in a cigar-shaped trap with lo^ = 2TT x 830 
Hz, iVy = 2tt X 415 Hz, and uiz — 2Tr x 22 Hz. Wc use the 
relation f « l.bT/Tp with kTp = {SNy/^tiLu to convert 
the effective temperature T obtained by fitting the den- 
sity profile to the real temperature [J, |9| . The calculated 
frequency and damping from ([5]) and ([9]) is plotted as 
solid black lines. The number of trapped atoms is taken 
to be = 4 X 10^. The dashed black lines in Fig. [T] are 
the hydrodynamic and coUisionless limits given by (llOp 
and (fTT|) respectively. We see that there is good overall 
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FIG. 1; (Color online) The frequency and damping of the 
scissors mode in the unitarity limit. The x 's are experimental 
values from Q. The black lines are the theoretical values 
obtained from the solutions to ([5]), while the dashed lines 
indicate the asymptotic results (|10|l and . The green dash- 
dotted curve is the damping obtained using (|12m3p and (|15p 
with u! — uih- 



agreement between theory and experiment both for the 
frequency and damping rate as a function of T. For low 
temperatures, the measured damping is somewhat larger 
than that obtained for our classical model. Since our 
model does not include Fermi blocking effects and super- 
fiuidity, this discrepancy is not surprising. Note that the 
theoretical curves contain no fitting parameters. 

Due to the anharmonicity of the trap, the observed 
frequency in the coUisionless limit is slightly reduced from 
PT|) [lo|, an effect which is not included in our model, 
which assumes the trap to be perfectly harmonic. 

We also plot in Fig. [T] as a dash-dotted green line the 
calculated damping rate obtained from hydrodynamics, 
using the approach described below. When the damping 
is small, the damping rate is determined by the ratio 
between the rate of loss of mechanical energy and the 
energy itself. For the velocity field given by ^ , the time 
average of the mechanical energy, -Emcch, is equal to 



mcch/ 



d^rmn{r)v^ (r) 



N . 



-kW 



1 



1 



(12) 

since in the classical limit the density is proportional to 
exp(— y(r)/fcr). The time average of the rate of loss of 
mechanical energy is (-Emcch) = ~2b^ J d^rrj, from which 
the damping rate is obtained as 



F = 



mcch / 



(13) 



2(-£'mcch) 

In the unitarity limit in the classical regime, a variational 
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calculation which is accurate to < 1% yields 



15 {mkTf/^ 
327^ P 



(14) 



for the viscosity. This is independent of density and the 
integration must be cut off close to the surface of the 
cloud where the gas is no longer hydrodynamic [l^ . A 
cut-off can also be introduced by using the real part of the 
(complex) dynamical viscosity 77 (w) = 77/[l — icL'r^(r)] [isj 
giving 



mccliy 



V 



l+c^V,(r)2 



where in the classical limit 



n{v)kT Nu) 



4.17 fkT\' 



.V(r)/kT 



(15) 



(16) 



As can be seen in Fig. [U this procedure leads to results 
which are qualitatively correct also in the coUisionless 
regime. It is interesting that the expression (|15p yields a 
damping rate which is exact in the coUisionless limit for 
a trap with cylindrical symmetry (co^ = ^y) provided lu 
is set equal to its value in the coUisionless limit, uj = ojci- 
For a general trap geometry, (fT5)) with lo = Wci yields a 
damping rate in the coUisionless limit which differs from 
(fTTj) by the factor Sa^{l + a)~2(l + a^y^ where a = 
LOy/uJx- For the trap parameters of the experiment in 
Ref. 0, the factor is « 0.71. 

We now examine the transition between hydrodynamic 
and coUisionless dynamics as a function of interaction 
strength. As in Ref. 01, we define the temperature 
where the damping is maximum as marking the transi- 
tion between the hydrodynamic and coUisionless regimes. 
In Fig. [21 we plot as a solid line Th calculated from ^ 
and ([9]) as a function of the interaction strength kpo,- 
The gas can be regarded as hydrodynamic below and col- 
lisionless above the line. We also plot the experimentally 
determined Th from Ref. As expected, the tran- 
sition between coUisionless and hydrodynamic behavior 
occurs at lower temperatures with decreasing interaction 
strength. There is good agreement between theory and 
experiment close to the unitarity limit. However, the cal- 
culations do not bring out the observed surprisingly steep 
decrease in Th with increasing values of — (fc_Fa)~^ on the 
BCS side of the resonance. This could be due to Fermi 
blocking effects making the system less hydrodynamic. 

The calculations reported here do not take into account 
Fermi blocking or superUuidity. In the weak-coupling 
limit (/cfIoI <C 1) it is straightforward to include Fermi 
blocking effects in the collision integral as described in 
Ref. la- The effect of superfluidity was discussed in 
Ref. [TJ for a uniform superfluid. The viscosity, which is 
associated with the motion of the normal component, de- 
creases below the transition temperature because of the 
energy gap in the spectrum of elementary excitations. 




FIG. 2: (Color online) The crossover temperature Th between 
hydrodynamic and coUisionless dynamics as a function of in- 
teraction strength. The x 's are experimental values from [3] 
and the soUd line is obtained as described in the text. 



The corresponding change in the viscous relaxation rate 
is numerically quite small (see Eq. (28) of Ref. Q). For 
a trapped superfluid the energy gap depends on position, 
and it is therefore much more difflcult to give a quanti- 
tative account of the effects of superfluidity on the mode 
frequencies. 



CONCLUSION 

We have analyzed the scissors mode of an interacting 
Fermi gas in the classical regime. By taking moments 
of the Boltzmann equation, we have calculated the fre- 
quency and damping of the mode both as a function of 
temperature and interaction strength. The calculation 
reproduces the main features of the recent experimental 
findings and can be used to identify non-classical effects 
for strongly interacting Fermi gases. 
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